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Abstract 

We study the sensitivity of wave travel times to steady and spatially homoge- 
neous horizontal flows added to a realistic simulation of the solar convection 
performed by Robert F. Stein, Ake Nordlund, Dali Georgobiani, and David 
Benson. Three commonly used definitions of travel times are compared. We 
show that the relationship between travel-time difference and flow amplitude 
exhibits a non-linearity depending on the travel distance, the travel-time def- 
inition considered, and the details of the time-distance analysis (in particular, 
the impact of the phase-speed filter width). For times measured using a Gabor 
wavelet fit, the travel-time differences become nonlinear in the flow strength for 
flows of about 300 m s _1 , and this non-linearity reaches almost 60% at 1200 m 
s _1 (relative difference between actual travel time and expected time for a linear 
behaviour). We show that for travel distances greater than about 17 Mm, the 
ray approximation predicts the sensitivity of travel-time shifts to uniform flows. 
For smaller distances, the ray approximation can be inaccurate by more than a 
factor of three. 
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1. Introduction 

Since the introduction of time-distance helioseismology (Duvall et at, 1993) this 
method has been applied to investigate various subsurface physical properties 
of the quiet Sun (e.g. Kosovichev and Duvall, 1997), especially supcrgranula- 
tion (e.g. Duvall and Gizon, 2000; Zhao and Kosovichev, 2003; Hirzberger et 
at, 2007), and sunspots (e.g. Kosovichev, Duvall, and Scherrer, 2000; Zhao, 
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Kosovichcv, and Duvall, 2001; Couvidat, Birch, and Kosovichev, 2006). In these 
studies, the main emphasis has been on mapping both the subsurface material 
flows and deviations of the sound speed from a reference solar model. Travel 
times of wavcpackets (whether from acoustic or surface gravity waves) are the 
basic measured quantities that are then inverted to produce estimates of the 
flow velocity and sound-speed perturbations. Travel times are measured from 
temporal cross-covariances. Over the years different definitions of these travel 
times have been employed {e.g., Kosovichev and Duvall, 1997; Gizon and Birch, 
2002; Gizon and Birch, 2004). For a given cross-covariance, these definitions, in 
general, yield different travel times. It is therefore important to understand how 
these definitions perform in simple situations where the subsurface conditions are 
known. In this paper we study the sensitivity of these definitions to horizontal, 
steady, and uniform flows. We use vertical velocity data from a realistic nu- 
merical simulation of solar convection (Benson, Stein, and Nordlund, 2006). We 
introduce steady and spatially uniform horizontal flows to this vertical velocity 
data cube and measure the perturbations in the travel times induced by these 
flows. We also compare the observed relationship between travel-time differences 
and flow velocities with the relationship predicted by ray approximation kernels 
(Kosovichev and Duvall, 1997) commonly used for the inversions. We limit our- 
selves to the study of acoustic modes (p modes), neglecting the surface-gravity 
modes (/ modes) . The /-mode case was studied in detail by Jackiewicz et al. 
(2007). In Section 2 we briefly present the numerical simulation that we used 
for this work. In Section 3 we remind the reader of the fundamentals of the 
time-distance formalism. In Section 4 we present our results, and we conclude 
in Section 5. 

2. Numerical Simulation 

We work with simulated linc-of-sight velocity data [<fi(r, t) where r is the horizon- 
tal position vector and t is time] obtained from a numerical three-dimensional 
simulation of convection in the upper solar convection zone. The numerical slab 
is 96 Mmx96 Mm wide and 20 Mm deep. The spatial resolution is 96 km 
horizontally and varies from 12 to 75 km vertically. The temporal resolution is 10 
seconds and we use 8.5 hours of simulated data. This simulation was performed 
by Robert Stein, Ake Nordlund, Dali Georgobiani and David Benson, and is 
publicly available on the website sha.stanford.edu. A very similar simulation (48 
Mmx48 Mm horizontally) has already been extensively analyzed by Georgobiani 
et al. (2007), Braun et al. (2007), and Zhao et al. (2007). The latter applied the 
time-distance formalism to these data to study, among other things, the validity 
of the inversion procedures. A detailed presentation of the code that produced 
these data can be found in Benson, Stein, and Nordlund (2006) and Stein and 
Nordlund (2000). Here we only work with the vertical velocity at the surface of 
the numerical simulation. The data cube we use is rebinned four times in the 
horizontal direction (yielding a final spatial resolution e?a;=384 km) and six times 
in time (final temporal resolution dt=60 seconds) . One advantage of using this 
numerical simulation in place of solar data is that horizontally the boundary 
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conditions are periodic. Therefore it is possible to simulate the impact of a 
steady horizontally uniform flow by simply shifting the vertical velocity maps in 
the horizontal direction at each time-step, using the shift theorem in the Fourier 
domain. Simulating an advection flow on actual solar data. e.g. on data from the 
SOI/MDI instrument (Scherrer et at, 1995), is possible (Jackiewicz et at, 2007) 
but more complicated and under normal conditions the simulated flow is non- 
uniform in latitude. The power spectrum of the vertical velocity <f>(r, t) from 
the numerical simulation is very similar to the actual solar oscillations power 
spectrum as determined by SOI/MDI (Georgiobani et al., 2007). Therefore we 
expect that the results we obtain here should be fairly representative of the 
actual solar case. 



3. Time-Distance Formalism: A Brief Overview 

3.1. Cross-Covariances 

We use the time-distance measurement procedure described in Couvidat, Birch, 
and Kosovichev (2006), except that we work with a center-quadrant geometry 
instead of a center-annulus one. We only deal with acoustic waves (p modes) 
and filter out the surface gravity waves (/modes) (see Jackiewicz et ai, 2007, 
for a recent study of the travel times of /modes). The measurement of the travel 
times begins with the computation of cross-covariances between two observation 
points ri and T2- To reduce the noise level of the cross-covariances, the data 
cubes are phase-speed filtered (Duvall et at, 1997). The data cube cj>(r,t) is 
filtered in the Fourier domain using a Gaussian filter F[k, w;A), where k is 
the horizontal wavenumber and u is the temporal angular frequency, for each 
travel distance A = ||r2 — ri||. An /-mode filter is also applied to remove the 
/-mode ridge. The point-to-point cross-covariances are averaged over quadrants 
(arc of 90°) of radius A centered on the point n. The quadrant geometry (two 
separate pairs of quadrants: East -West and North -South) is used to retrieve 
the subsurface horizontal flow velocity. Such point-to-quadrant cross-covariances 
are computed for four central distances (A) ranging from 6.2 Mm to 16.95 Mm. 
We do not select larger A because there is a small amount of wave reflection at 
the bottom boundary of the numerical simulation (Zhao et a/., 2007) that seems 
to alter the shape of the cross-covariances. The larger the A the deeper the 
wavepackets reach, and the more they are affected by these reflected waves. Ray 
theory shows that for A = 16.95 Mm, the lower turning point of the wavepackets 
is located at a depth of about 4.5 Mm, far from the bottom boundary. To further 
increase the signal-to-noise ratio of the cross-covariances, we average them in 
groups of five distances (A) for each central distance. For instance, the final 
cross-covariances at the central distance A = 6.2 Mm are averages of five cross- 
covariances computed at distances ranging from A = 3.7 to 8.7 Mm (see details 
in Couvidat, Birch, and Kosovichev, 2006). The values of A and properties of 
the phase-speed filters used in this paper are summarized in Table 1. 
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Table 1. Parameters for the time-distance analysis: distances A = ||r2 — ri||, range 
of A, properties of the Gaussian phase-speed niters applied to the Fourier transform 
of the vertical-velocity data cube (dq is the central phase speed and FWHM is the full 
width at half maximum), centers of the time window used to measure the travel times, 
and selected phase travel times (t p ) of the reference cross-covariances 



A (Mm) 


range (Mm) 


vo (km s 1 ) 


FWHM (km s" 1 ) 


to (min) 


t p (min) 


6.2 


3.7-8.7 


12.77 


6.18 


19.0 


10.48 


8.7 


6.2-11.2 


14.87 


6.18 


21.4 


12.85 


11.6 


8.7-14.5 


17.49 


6.18 


24.4 


19.65 


16.95 


14.5-19.4 


24.82 


9.09 


28.7 


29.05 



3.2. Travel-Time Definitions 

The point-to-quadrant cross-covariances are used to determine the travel times 
of the acoustic wavepackcts. Traditionally, cross-covariances have been fit with 
a Gabor wavelet defined as (Kosovichev and Duvall, 1997): 

(c 2 \ 
~r (t-Tg) 2 ) (!) 

where A is the amplitude of the wavelet, too is the central temporal frequency 
of the wavepacket, 5uj is a measure of its frequency width, r p is the phase travel 
time, and r g is the group travel time. It is customary to use only phase travel 
times with the time-distance formalism, because group travel times are more 
difficult to retrieve accurately (the center of the Gaussian envelope is usually 
poorly estimated) and therefore group travel-time maps are significantly noisier 
than phase travel-time maps. 

Recently, two more definitions of the wavepacket travel times based on studies 
in geophysics were added. Gizon and Birch (2002) [hereafter GB02] define the 
travel times r±(r, A) as: 

t± (r, A) = arg min t X± (r, A, i) (2) 
where X±(r, A, t) are the following functions: 

X ± (v, A, t)=dtJ2 f(±t')[C(r, A, t') - C rcf (A, t' T t)] 2 (3) 
v 

where f(t) is a one-sided window function that selects the first-bounce skip of 
the cross-covariance (the skip corresponding to waves travelling directly from 
the center to the quadrant) on the positive-time or negative-time branches of 
the cross-covariance. The same window function is also used when fitting the 
cross-covariance by a Gabor wavelet even though it does not appear explicitly in 
Equation (1). The function C rof (A,t) is the reference cross-covariance. Usually 
C (A, t) is either obtained from a model of the solar oscillation power spectrum, 
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or is a spatial average of the measured cross-covariances over a region of the quiet 
Sun. 

Lastly, Gizon and Birch (2004) [hereafter GB04] proposed another definition 
of the travel time that is linear with respect to the cross-covariance: 

r± (r, A) = dtJ2 W ± (A, t) [C(v, A, t) - C rcf (r, A, t)] (4) 
t 

where the W± are weight functions: 

W ± (A,t)= (5 ) 

and C (A,t) denotes the time derivative of the reference cross-covariance. 

As emphasized by the use of the ± symbol, for each of the three definitions 
of travel times we can measure r on the positive-time branch (for wavepack- 
ets travelling in a specific direction, e.g. from the South to the North) and 
on the negative-time branch (for wavepackets travelling in the opposite direc- 
tion). Therefore we measure four different kinds of travel time: r\v(r, A) is 
the travel time of wavepackets traveling westward over a horizontal distance 
A centered on r; t e (r, A) is for eastward traveling waves; rs(r, A) is for south- 
ward traveling waves; and finally tn(i", A) is for northward traveling waves. The 
corresponding travel-time perturbations are produced by subtracting the travel 
time obtained from the reference cross-covariance. For instance, <5tn(i", A) = 
tn(i", A) — TN ;rc f(A) where TN ;ro f(A) is the travel time of the reference cross- 
covariance. Here, C* rof (A, t) is obtained by averaging over r the center-to-annulus 
cross-covariance. In this paper we focus on 5tns(i", A) = <5tn(i", A) — <5rs(r, A), 
the North- South difference travel-time perturbation, that is, in a first approx- 
imation, only sensitive to horizontal flows in the North -South direction. With 
this convention a northward flow will result in a negative <5tns- 



4. Results 

We added southward flows to the velocity data cube from the numerical simula- 
tion by applying the shift theorem in the Fourier domain. At each time step (t) 
we shift the velocity maps <p(r, t) in the y direction (from North to South) by 
the amount Sy = ut where u is the amplitude of the imposed flow. Again, the 
advantage of using the simulation for this purpose is that its boundary conditions 
are periodic. We tested 12 flow amplitudes: u = 0, 100, 200, 300, 400, 500, 600, 
700, 800, 900, 1000, and 1200 m s" 1 . Even in the absence of any additional flow 
(case u = m s _1 ) the simulation already includes a weak northward flow. This 
flow is horizontally non-uniform and its velocity varies slightly with depth, with 
a maximum amplitude of the order of 35 m s . The simulation also includes 
a horizontally non-uniform eastward flow that is somewhat stronger (with an 
average velocity of about 150 m s _1 in the first few Mm underneath the surface) 
and is strongly depth dependent. The presence of this stronger eastward flow 
explains why we favor the North - South direction in the results presented here. 
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However we also tested the addition of eastward flows and their impact on the 
East -West travel-time difference, and obtained similar results. In the rest of 
this paper we use two kinds of Gaussian phase-speed filters when computing 
the cross-covariances: those whose parameters are listed in Table 1 (hereafter 
standard phase-speed filters), and those whose FWHM was multiplied by four 
(hereafter broad phase-speed filters). 

4.1. Impact of the Horizontal Flows on the Cross-Covariances 

The impact of the additional uniform horizontal flows on the spatially averaged 
cross-covariances is shown in Figure 1 (at A = 8.7 Mm). Because the flows 
are southward, the spatially averaged cross-covariance for the northward-going 
wavepackets (hereafter northward cross-covariance) is shifted toward larger times 
compared to the reference cross-covariance; the northward-going wavepackets are 
slowed down by the southward flows. The reverse is true for the southward-going 
wavepackets. However flows do not result only in a shift of the cross-covariances. 
Indeed, when we fit the cross-covariances by Gabor wavelets (Equation (1)), 
other differences unfold. 

i) The central frequencies (u>o) of the wavepackets are also different for the 
northward and southward cross-covariances (Figure 2). For instance, at A = 8.7 
Mm, in the presence of a 400 m s _1 southward flow, loq = 3.936 x 2tt 10 -3 rad 
s _1 for the spatially-averaged northward cross-covariance, while a>o = 4.247 x 2-7T 
10 -3 rad s _1 for the southward cross-covariance. The change in uiq is generally 
consistent with a Dopplcr shift Suiq described by Suiq = k • u and resulting from 
a flow u. However at short distances (A = 6.2 and 8.7 Mm), and for flows faster 
than about u = 600 m s _1 , 5u>q becomes non-linear in u. This non-linearity 
is probably mainly the result of the application of Gaussian phase-speed filters 
to the Fourier transform of the data cube. The k • u Dopplcr shift changes 
the location of the acoustic-mode ridges in regard to the filters and affects the 
power distribution inside these filters depending on the ridge distance to the 
filter center. Moreover, this change in the power distribution is different for 
the positive and negative (in term of ui) branches of the acoustic-mode ridges, 
explaining why the luq of northward and southward cross-covariances are not 
symmetrical around the ujq value at rest (u = 0). With the broad phase-speed 
filters at A = 6.2 Mm, the acoustic wavepackets have a lower central frequency 
(wo) and a lower central wavenumber (fc) where there is more acoustic power 
than with the standard filters. Therefore, with the broader filters, 6u>o is smaller 
for a given velocity u, because the average k of the wavepacket is smaller. Finally, 
5uo becomes non-linear in u earlier with the broad filters than with the standard 
ones, at least at short A: this might be because the average k is more dependent 
on the flow amplitude than in the case of the standard filters. 

ii) The frequency width of the wavepacket varies (Figure 3): for u = 400 
m s -1 , Suj = 0.74 x 2ir 10~ 3 rad s _1 for the northward cross-covariance and 
Su> = 0.91 x 2-7T 10 -3 rad s _1 for the southward one. A change in the wavepacket 
frequency width is expected because the Doppler shift Su>q depends on k, but 
as with the change in wq, the interplay between the w-shifted acoustic ridges of 
the power spectrum and the Gaussian phase-speed filters probably explains the 
details of this frequency width change. 
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Hi) The maximum amplitude of the cross-covariances is also affected by the 
presence of a flow, and again this can partly be explained by a change in the 
amount of acoustic power selected by the phase-speed filters due to the Doppler 
shift. 

All of these changes affect the northward and southward cross-covariances 
differently This can also been seen in Figure 4 which shows the cross-covariances 
as a function of both the time and the flow amplitude for standard and broad 
phase-speed filters. 

There is a clear impact of the FWHM of the filter on the shape of the cross- 
covariances (and, consequently, on the values of the travel times returned by the 
definitions studied here) . This impact has been mentioned regarding some Gabor 
wavelet parameters (Figures 2 and 3). Another consequence of the phase-speed 
filters, not modeled by the Gabor wavelet, is the presence of a filter artifact (e.g. 
Couvidat et al, 2005). Indeed, the cross-covariances for the distances studied in 
this paper have three main components: a filter artifact (at small times \t\), a 
first-bounce skip (waves travelling directly between the two observations points), 
and a second-bounce skip (waves being reflected once at the surface in between 
ri and r2). At A = 6.2 and 8.7 Mm, the filter artifact, the first-bounce skip, and 
the second-bounce skips are entangled. It appears in Figures 1 and 4 that the 
filter artifact is the part of the cross-covariance that is the least affected by the 
presence of a flow: the location and amplitude of its peaks do not vary much. 
Standard phase-speed filters engender filter artifacts whose amplitude is larger 
than those arising from broad filters. Fast, southward flows push the first-bounce 
skip of the southward cross-covariances toward lower times \t\. Therefore, the 
peaks of this first-bounce skip become more and more altered by the presence of 
the filter artifact. On the other hand, the use of broad phase-speed filters reduce 
the amplitude of the artifact and the peaks of the first-bounce skip are less 
distorted. An example of the time window [/(£)] used to select the first-bounce 
skip is shown in Figure 1. For the travel-time measurements this time- window 
is always centered at the same time (independent of the flow amplitude) . 

Figure 5 shows the power-spectral density, and the real and imaginary parts 
of the Fourier transform of the northward and southward spatially averaged 
cross-covariances in the presence of a 1000 m s _1 southward flow, for the stan- 
dard and broad phase-speed filters at A = 8.7 Mm. The differences mentioned 
earlier in central frequency, frequency width, and total power of the wavepacket 
between northward and southward cross-covariances are clearly visible. It con- 
firms that the asymmetry between southward and northward cross-covariances 
is mainly a result of the interaction between the phase-speed filters and the 
wavepackets. Overall, broader phase-speed filters seem preferable when com- 
puting the cross-covariances, since narrow filters introduce more asymmetries 
between northward and southward cross-covariances and significantly alter the 
shape of the first-bounce ridge at short distances A. This could lead to spurious 
interpretations regarding physical phenomena like subsurface flows or sound- 
speed velocity perturbations, when the resulting travel times are inverted with 
sensitivity kernels that do not take into account the details of the travel-time 
measurement technique. 
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4.2. Definition of the Phase Travel Time With the Gabor Wavelet 

As mentioned in Section 3.2. when fitting a temporal cross-covariance with a 
Gabor wavelet we obtain two kinds of travel times: a group r g and a phase r p 
travel time. A well-known issue with r p is that it is not uniquely defined {e.g. 
D'Silva, 2001). Indeed any time congruent to r p modulo 2-k/loq is an equally 
valid phase travel time because the Gabor wavelet is the same for any of these 
numbers. 

In return, the lack of uniqueness of t p means that the value of (5tns(i",A) 
engendered by a southward flow is also not unique and depends on the value of 
the selected reference phase travel time. Suppose that we select 12.85 minutes 
as the phase travel time of C ref (A,t) for A = 8.7 Mm, and we add a 200 m 
s _1 southward flow to the numerical simulation data cube. Then the Gabor 
wavelet fitting code returns tn = 12.917 minutes and r§ = 12.781 minutes, 
resulting in <5tns = tns = 8.15 seconds. If we select as reference phase travel 
time 12.85 + 2tt/ujq = 16.95 minutes instead, the fitting code returns tn = 
12.917 + 2-7r/cj ;N = 17.074 minutes and r s = 12.781 + 2tt/uj . s = 16.794 minutes, 
resulting in (Stns = 16.79 seconds. Indeed, the too parameter for the northward 
cross-covariance is luq.^ — 4.009 x 2ir 10~ 3 rad s _1 while for the southward 
cross-covariance it is ujq.^ = 4.153 x 2ir 10~ 3 rad s _1 . This difference in luq is a 
direct consequence of the 200 m s^ 1 flow. This example shows that the same flow 
velocity can produce different values for <5rNs( r , A) (here 8.15 and 16.79 seconds). 
This is a major difference between the three travel-time definitions studied here: 
GB02 and GB04 return a uniquely defined travel time, unlike the phase travel 
time of the Gabor wavelet. Figure 6 shows the impact of the reference phase 
travel time on i5tns as a function of the southward flow velocity, at A = 16.95 
Mm. 

Therefore, when using the Gabor wavelet definition it is necessary to define a 
criterion regarding which values of the phase travel time are to be selected. It is 
customary (e.g. J. Zhao, private communication) to choose as phase travel times 
those whose value is the closest to — and if two phase travel times are equally 
distant, smaller than — the corresponding group travel times (which is uniquely 
defined). It can be argued that selecting r p close to r g makes sense because 
that is where the cross-covariance peaks have the largest amplitude and are the 
furthest away from both the filter artifact and the second-bounce skip, therefore 
making the determination of their location more accurate. At A = 8.7 Mm the 
reference group travel time is about equally distant from 12.85 and 16.95 minutes. 
Therefore the phase-travel time that we select is 12.85 minutes. Trying to fit first- 
bounce skip peaks located at shorter times is a difficult task as can be seen in 
Figure 4. In this figure, the dashed lines representing the phase times do not 
follow the peaks of the cross-covariances for times |t| < 10 minutes, confirming 
that these peaks are mainly filter artifacts. This criterion also has the advantage 
that the selected t p is close to the phase travel time t p = J T k t /uds (T is the 
ray path, k t is the total — not just horizontal — wavenumber, and ds is the 
element of path length) returned by the ray-path approximation at short A 
distances (11.3 minutes at A = 8.7 Mm, when applied to the solar model S of 
Christensen-Dalsgaard et al, 1996). At large A this is not the case anymore: for 
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instance at A = 16.95 Mm, r g = 26.10 minutes and t p = 16.51 minutes according 
to the ray-path approximation based on model S, while the equal-to-or-smaller- 
than criterion implies a reference r p = 23.51 minutes. The reference phase travel 
times we select for the rest of this paper when fitting a Gabor wavelet are listed 
in the last column of Table 1 . 

4.3. Comparison of Definitions of Travel Times 

Figure 7 shows the values of the North- South difference travel-time pertur- 
bation (5tns as a function of the amplitude of the horizontal uniform flow. All 
three travel-time definitions arc shown. As mentioned in the previous section, the 
Jtns for the Gabor wavelet fit can be changed by selecting a different reference 
t p , and thus make them closer to the ones obtained with GB02. The 5tns(A) 
displayed are the expected values of Gaussian functions fitted to the histograms 
of <5TNs( r ;A). Fitting a Gaussian seems more appropriate than reporting the 
spatial average of <5tns because of the presence of numerous misfits with the 
Gabor wavelet and GB02 definitions (the larger the flow, the more misfits the 
codes return), and the distribution of 5tns(f, A) at a given A is very close to a 
Gaussian distribution. The presence of an intrinsic horizontal flow in the numer- 
ical simulation explains why none of the <5tns (A) is equal to zero in the absence 
of any additional flow. Both the values of <5tns measured from cross-covariances 
obtained with the standard and broad phase-speed filters are shown. 

There is a strong dependence of (*>tns on the width of the phase-speed filter at 
small distances A. This is hardly surprising considering the significant impact of 
the phase-speed filters on the shape of the cross-covariances. Overall, the three 
travel-time definitions give significantly different answers to the question of what 
the travel time of a wavepacket is. At low flow velocities, GB02 and GB04 agree, 
as expected since GB04 can be seen as a linearization of GB02. However a non- 
linearity with respect to the flow velocity sets in early for GB04 (at about 200 m 
s _1 at A = 6.2 Mm) and the discrepancy with GB02 becomes increasingly larger 
for larger flows. Indeed, GB04 appears to be inadequate for the detection of flows 
larger than about 700 m s _1 (at A = 6.2 Mm and with standard phase-speed 
filters). We note, however, that GB04 was designed to be applied to quiet-Sun 
data only. Even though at short A all three travel-time definitions exhibit a non- 
linearity for the flows tested here, GB04 is the definition departing the earliest 
from linearity For the three definitions the non-linearity is larger with narrower 
filters (except at A = 16.95 Mm with GB04). This emphasizes the role of the 
phase-speed filter in the creation of this non-linearity. At larger A, GB02 and 
the Gabor wavelet fit return similar values for the travel times. 

With the standard phase-speed filters, no travel-time definition returns a 
unique relation between <5tns and the flow velocity at A = 6.2 Mm within the 
range of flow amplitudes tested. Indeed, a measured <5tns = 17 seconds might 
be interpreted as arising from a 1200 m s _1 or a 500 m s _1 flow according 
to the Gabor wavelet fit. This means that the inversion procedure will not 
be able to detect flows of the order of 1200 m s _1 , and will instead assign a 
<5t"ns = 17 seconds value to a 500 m s _1 flow, thus underestimating the true 
flow velocity. Therefore flows faster than s» 1000 m s _1 in the shallow layers 
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underneath the solar surface are impossible to retrieve if we are not careful 
with the details of the time-distance analysis (width of the phase-speed filters). 
However such high flow velocities are essentially never present in solar data 
except in untracked data cubes (cubes not corrected for the solar rotation). 
Figure 8 shows the derivatives of <5tns with respect to the flow amplitude. To 
display smooth derivatives, the 5tns were first fitted with a polynomial of order 
three and the derivatives were computed from these polynomials. At A = 6.2 Mm 
the relation between Jtns and flow velocity becomes not linear for flows faster 
than roughly 300 m s _1 for the Gabor wavelet and GB02 definitions and with 
standard phase-speed filters. Broader filters delay the onset of non-linearity to 
about 600 - 700 m s" 1 . However, at A = 11.6 and 16.95 Mm, GB02 and the 
Gabor wavelet are approximately linear for the range of flow velocity tested, 
except maybe at 1200 m s _1 . GB04 is never linear in flow velocity over the 
entire range tested in our study. At A = 6.2 Mm, and for standard phase-speed 
filters, the derivatives of <5tns become negative for a ratio of the flow velocity 
(v) to the phase-speed filter central phase speed (vq) larger than 6.5%. At this 
point the definitions of travel times become inadequate because they cannot 
discriminate between fast and slow flows, as previously mentioned. GB04 is the 
definition for which the derivatives turn negative at the slowest flows, followed 
by the Gabor wavelet, and GB02. The ratio v/vo at which the derivative of 5tns 
becomes negative becomes smaller and smaller with increasing A for GB04. 

The non-linearity and disagreement among travel-time definitions seem to 
decrease as A increases, which is the opposite of what occurs in Figures 3 and 5 
of Jackiewicz et al. (2007), for a travel-time difference based on GB04. However, 
Jackiewicz et al. (2007) used ridge filtering (f-mode filter) and surface-gravity 
waves, while here we use phase-speed filtering and acoustic waves. Moreover 
because of differences in the physical models, the oscillation power spectra are 
also different. Therefore, it is difficult to draw any firm conclusion from this 
possible discrepancy between the two papers. 

Lastly, we remind the reader that only GB04 is linear with respect to the cross- 
covariance. Therefore, the travel time of the spatially-averaged cross-covariance 
is not equal to the spatial average of the travel times with GB02 and the Gabor 
wavelet. We find that the difference between the average of the travel times and 
the travel time of the average cross-covariance increases with the flow velocity. 
This is in part due to the increasing number of misfits we encounter with these 
two travel-time definitions when the flow velocity increases. 

4.4. Frequency Dependence of the Travel-Time Differences 

Following Braun and Birch (2006) and Couvidat and Rajaguru (2007) we applied 
Gaussian frequency filters to the Fourier transform of the data cube, on top of the 
phase-speed filters, to select three specific frequency bands. These three bands 
are centered onv = 3.5, 4.0, and 4.5 mHz. The FWHM of these Gaussian filters 
is 1 mHz. Figure 9 shows the frequency dependence of <5tns as a function of flow 
velocity for cross-covariances obtained with the standard phase-speed filters. 
Figure 10 shows the same but for broad phase-speed filters. Again the impact 
of the phase-speed filter width is clearly visible in these plots. The travel-time 
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difference perturbation exhibits a strong dependence on the frequency. For the 
Gabor wavelet and GB02 the higher central frequency of the wavepacket leads 
to larger <5tns- 

We expect to observe a frequency dependence of the travel times: even though 
the argument that wavepackets with the same central phase-speed (w/fc) follow 
approximately the same ray path is used to justify the application of phase-speed 
filters prior to the computation of cross-covariances, the ray theory shows that 
in actuality there can be a bit of a difference in the trajectories. Thus, both the 
upper and lower turning points of a wavepacket change with its central frequency. 
For instance at A = 16.95 Mm, with the solar model S, the upper turning point is 
located about 4.8 km deeper for a wavepacket of central frequency v — 3.5 mHz 
than for a wavepacket of central frequency v = 4.5 mHz. The lower turning point 
changes by more. However the reversal in the sign of <5tns observed for fast flows 
at A = 6.2 Mm and a frequency of 3.5 mHz cannot be explained by a small 
change in the ray-path of the wavepacket. The complex interplay between the 
phase-speed filter and the wavepacket seems to be the main reason behind this 
strong non-linearity at low frequency. This result confirms again the interest in 
broad enough phase-speed filters in time-distance analysis, and that we cannot 
ignore the details of the travel-time measurement technique when interpreting 
the travel times. 

4.5. Comparison with the Ray-Path Approximation 

To invert the travel times and retrieve the physical perturbations to which these 
times are sensitive, sensitivity kernels (i.e. a model for the linear dependence 
of travel-time shifts on perturbations to a solar model) must be employed. Two 
main approximations have been used to produce these kernels: the ray-path 
(Kosovichev and Duvall, 1997) and Born approximations (Birch et al, 2001; 
Gizon and Birch, 2002; Birch and Gizon, 2007, in the context of time-distance 
helioseismology) . The ray-path approximation (or geometric optics approxima- 
tion) is valid in the limit of small wavelength while the Born approximation takes 
into account finite wavelength effects. In the ray approximation, the travel-time 
difference St caused by a flow v is (Kosovichev and Duvall, 1997): 



where n is the unit vector tangent to the ray path, and c is the sound speed. 
In both the ray and Born approximations, the relationship between physical 
perturbation and travel-time perturbation is assumed to be linear. As can be 
seen in Figure 7, the assumption of linearity with the amplitude of the flow 
is a major issue when the non-linearity of the travel times sets in. This shows 
how far the inverted flow velocities could be from the actual flows, for flows 
with amplitudes above a few hundred m s _1 . The thick black lines on all four 
panels show the variation of <5tns(A) — as a function of the flow velocity - 
predicted by Equation (6) applied to model S (which is close to the background 
model of the numerical simulation). At A = 6.2 Mm and for standard phase- 
speed filters, the ray-path kernels seem to overestimate by a factor larger than 




(6) 
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three the sensitivity of travel times to flows. None of the travel-time definition 
returns results close to the ones predicted by this ray-path approximation. A 
way to improve the agreement between predicted and measured travel times 
with the Gabor wavelet could be to select larger reference phase travel times. As 
mentioned in Section 3.1, our cross-covariances are always averaged over a A- 
distance range to improve their signal-to-noise ratio. On the other hand, the ray- 
path kernels we compute here are for the specific A listed in Table 1 (no average 
over a range of A). Therefore we do not expect perfect agreement between the 
£tns(A) predicted by the ray-path kernels and the measured <5tns(A). However, 
we also computed cross-covariances that are not averaged over A for comparison, 
and the travel times that we measured on these cross-covariances are very close 
to the travel times measured on the averaged cross-covariances (the maximum 
difference between the two (5tns(A) occurred at flows faster than 600 m s _1 and 
was only about five seconds). 

5. Conclusion 

We added spatially-uniform steady southward flows of amplitudes ranging from 
to 1200 m s _1 to the vertical velocity data cube of a numerical simulation of solar 
convection. These southward flows produced a perturbation in the North - South 
difference travel times of acoustic wavepackets. We measured these travel times 
using three different definitions: the Gabor wavelet (Duvall and Kosovichev, 
1997), Gizon and Birch (2002), and Gizon and Birch (2004). The customary use 
of the phase travel time returned by the Gabor wavelet fit instead of the group 
travel time creates an ambiguity regarding the value of the difference travel-time 
perturbation engendered by a given flow velocity. This is a direct consequence 
of the lack of uniqueness of the phase travel time together with the fact that 
horizontal flows shift the central temporal frequency of the wavepackets. This 
ambiguity in the travel-time difference is a major issue especially when we try to 
interpret this difference using sensitivity kernels. Conversely, GB02 and GB04 
return uniquely defined travel times. At short distances A between observation 
points, all three travel-time definitions exhibit a non-linearity with respect to 
the flow amplitude. At these short distances, the details in the calculation of the 
cross-covariances, and especially the width of the phase-speed filters, prove cru- 
cial. Broader phase-speed filters may be preferable because they push the onset 
of non-linearity to larger flow amplitudes. GB04 is the definition of travel times 
that departs the earliest from linearity: the onset seems to be at flow amplitudes 
as low as about 200 m s" 1 . Therefore the use of GB04 can be problematic when 
not applied strictly to the quiet Sun. In addition, we note that supergranules 
have flow velocities that can exceed 200 m s _1 , and therefore the use of inversion 
procedures on travel-time differences obtained by the GB04 definition could 
return incorrect velocities (in particular underestimated velocities). GB02 and 
the Gabor wavelet also exhibit a sudden drop in the travel-time difference at 
flows faster than 1 km s _1 at A = 6.2 Mm with standard phase-speed filters. 
This means that fast flows near the solar surface will produce the same travel- 
time difference as slower flows. Flows of these amplitude are not common in 
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tracked data cubes where the effect of solar rotation has been largely removed. 
The non-linearity of the travel-time shifts is reduced for the case of broad phase- 
speed filters. At distances (A) larger than or equal to 16.95 Mm both GB02 
and the Gabor wavelet definitions produce travel-time differences that vary 
essentially linearly with respect to flow amplitude, within the velocity range 
studied here. If a non-linearity sets in. it is for flows much faster than 1.2 km 
s _1 . At distances less than about 17 Mm, the sensitivity kernels produced by the 
ray-path approximation fail to explain the sensitivity of travel-time differences 
to flows. 
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Figure 1. Left column: spatially-averaged northward (blue) and southward (red) cross-co- 
variances for different southward flow velocities, at distance A = 8.7 Mm, and for standard 
phase-speed filters. The flow amplitudes are, from upper to lower panels: 0, 200, 400, 600, 
800, 1000, and 1200 m s _1 . In black is the reference cross-covariance C (A, t). The temporal 
window applied when measuring the wavepacket travel times is also shown. Right column: the 
difference between the northward (blue) or southward (red) cross-covariance and C* ref (A,t). 
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Figure 2. Gabor wavelet parameter uiq (wavepacket central frequency) as a function of the 
southward flow velocity. Black lines are for the northward cross-covariances, red lines are for 
the southward ones. Solid lines are for the standard phase-speed filters, dashed lines are for 
the broad phase-speed filters. The measurement points are shown as diamonds. 
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Figure 3. Gabor wavelet parameter 8w (wavepacket width) as a function of the southward 
flow velocity. Black lines are for the northward cross-covariances, red lines are for the south- 
ward ones. Solid lines are for the standard phase-speed filters, dashed lines are for the broad 
phase-speed filters. The measurement points are shown as diamonds. 
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Figure 4. Cross-covariances as a function of time and flow amplitude, for a distance A = 8.7 
Mm, and for standard phase-speed filters (upper panel) and broad phase-speed filters (lower 
panel). Superimposed in black lines are the phase travel times r p of the Gabor wavelet fit 
selected for this paper (solid lines) and some other phase travel times r p ± k X 2tt/loq where k 
is an integer (dashed lines). Notice that the filter artifacts at roughly |t| < 10 min are largely 
independent of flow amplitude. 
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Figure 5. Power spectral density (left panels), real (middle panels), and imaginary (right 
panels) parts of the Fourier transform of the northward (solid lines) and southward (dashed 
lines) spatially averaged cross-covariances at A = 8.7 Mm in the presence of a 1000 m s~~ 1 
southward flow. Upper panels: standard phase-speed filter; lower panels: broad phase-speed 
filter. 
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Figure 6. Impact of the choice of the Gabor wavelet phase travel time on the relation be- 
tween North — South difference travel-time perturbations and the southward flow velocity, at 
a distance A = 16.95 Mm. Black lines are for the standard phase-speed filter and red lines 
are for the broad filter. Solid lines are for a reference r p = 29.05 minutes (standard filter) or 
29.25 minutes (broad filter), dashed line is for r p = 23.5 minutes, and dash-dotted lines are 
for t p = 17.96 minutes (standard filter) or 18.63 minutes (broad filter). The group travel time 
is T g = 27.69 minutes (standard filter) or 27.71 minutes (broad filter). 
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Figure 7. North — South difference travel-time perturbations as a function of the southward 
flow velocity. The solid lines are for the travel times measured on cross-covariances obtained 
with the standard width of the phase-speed filters, while the dashed lines are for those obtained 
with four-times wider phase-speed filters. The black lines are the results of the Gabor wavelet 
fit, the green lines are for Gizon and Birch (2002), and the red lines are for Gizon and Birch 
(2004). The perfectly linear curves are the results of the ray-path approximation. 
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Figure 8. Derivatives of the North -South difference travel-time perturbations with respect 
to the southward flow velocity. The solid lines are for the travel times measured on cross-co- 
variances obtained with the standard width of the phase-speed filters, while the dashed lines 
are for those obtained with four-times wider phase-speed filters. The black lines are the results 
of the Gabor- wavelet fit, the green lines are for Gizon and Birch (2002), and the red lines 
are for Gizon and Birch (2004). The horizontal solid lines show the results of the ray-path 
approximation. 
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Figure 9. North — South difference travel-time perturbations as a function of the southward 
flow velocity and when Gaussian frequency filters have been applied on top of the phase-speed 
filters. All the travel times were obtained with the standard width of the phase-speed filters. 
The solid lines are for a frequency filter centered on 3.5 mHz, the dashed lines are for 4 mHz, 
and the dot-dashed lines are for 4.5 mHz. The black lines are the results of the Gabor wavelet 
fit, the green lines are for Gizon and Birch (2002), and the red lines are for Gizon and Birch 
(2004). 
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Figure 10. North- South difference travel-time perturbations as a function of the southward 
flow velocity and when Gaussian frequency filters have been applied on top of the phase-speed 
filters. All of the travel times were obtained with phase-speed filters four times wider than the 
standard width. The solid lines are for a frequency filter centered on 3.5 mHz, the dashed lines 
are for 4 mHz, and the dot-dashed lines are for 4.5 mHz. The black lines are the results of 
the Gabor wavelet fit, the green lines are for Gizon and Birch (2002), and the red lines are for 
Gizon and Birch (2004). 
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